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Abstract 

We employ holographic techniques to study quantum quenches at finite temperature, 
where the quenches involve varying the coupling of the boundary theory to a relevant 
operator with an arbitrary conformal dimension 2 < A < 4. The evolution of the 
system is studied by evaluating the expectation value of the quenched operator and 
the stress tensor throughout the process. The time dependence of the new coupling is 
characterized by a fixed timescale and the response of the observables depends on the 
ratio of the this timescale to the initial temperature. The observables exhibit universal 
scaling behaviours when the transitions are either fast or slow, i.e., when this ratio is 
very small or very large. The scaling exponents are smooth functions of the operator 
dimension. We find that in fast quenches, the relaxation time is set by the thermal 
timescale regardless of the operator dimension or the precise quenching rate. 
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1 Introduction 

Recent advances in cold atom experiments have stimulated a vigorous research pro- 
gram into quantum quenches, processes in which the physical couplings of a quantum 
system are abruptly changed [1]. The basic motivation is to understand the organiz- 
ing principles governing the far-from-equilibrium behaviour of such systems. Although 
such quenches are well understood in the context of quantum mechanics [2], much less 
is known about such processes in quantum field theories. Theoretical progress has 
been made for a variety of systems, including two-dimensional conformal field theories, 
(nearly) free field theories and integrable models - e.g., see [3,4]. However, broadly ap- 
plicable theoretical techniques, which provide an efficient description of these quenches, 
remain to be found. 

Gauge/gravity duality [5] provides a remarkable new approach to studying certain 
strongly coupled field theories. Of course, in this framework, questions about the field 
theory are recast into questions about gravity in one higher dimension. These holo- 
graphic models seem to be especially well suited for the study of quantum quenches 
since, with relatively modest efforts, one is able to study strongly coupled quantum 
field theories, real-time processes and systems at finite temperature, as well as allow- 
ing for analysis in general spacetime dimensions. Hence holographic techniques have 
recently been applied to the study of quantum quenches [6] and the related issue of 
'thermalization' [7,8]. However, given the complexities of the bulk description of rapid 
changes in the boundary theory, numerical relativity is increasingly being applied to 
study these far-from-equilibrium processes [9,10] — see also [11]. 

In this paper, we extend the calculations presented in [12], in which the gauge/gravity 
duality was used to study 'thermal quenches' in a plasma of the strongly coupled 
jV = 2* gauge theory. More specifically, [12] studied the response of an initial ther- 
mal equilibrium state to variations of the coupling of the boundary theory to either a 
dimension two or three operator. The analysis was restricted to a high temperature 
regime where the calculations were carried out to leading order in m/T <^ 1. Here 
m is the relevant mass scale introduced by the new coupling. In the present case, 
we extend these holographic calculations to consider quenches made by coupling to 
a relevant operator with an arbitrary conformal dimension in the range 2 < A < 4. 
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The behaviour of the strongly coupled boundary theory in the present quenches is very 
similar to that found in [12]. In fact, for many of the results, we are able to identify 
a simple function of the conformal dimension which interpolates between the different 
cases which are explicitly studied both here and in [12]. For example, we find that in 
fast quenches, the increase in the energy density scales like {T^/ Aty^~'^, where Tj is 
the initial temperature and At is the timescale over which the new coupling is turned 
on — fast quenches are then those for which Ti/At ^ 1 . 

The remainder of the paper is organized as follows: In section 2, we describe the 
holographic model which is used to study our quenches and derive the gravitational 
equations that are to be solved. Next we examine solutions of these equations in 
section 3. In particular, by restricting our attention to the high temperature regime, 
we show that to leading order we only need to solve the linearized equation for the 
bulk scalar. We also consider the asymptotic boundary expansion for these solutions in 
Eddington-Finkelstein coordinates. In section 4, we translate the latter expansion to 
Feffer man- Graham coordinates, which are more suitable to stTidy physical observables 
in the boundary theory. In section 5, after finding the counterterms that renormalize 
the bulk action, we find expressions for the expectation value of the stress-energy tensor 
and the quenched operator in terms of gravitational variables. We also show that these 
observables obey the expected Ward identities. In section 6, wc identify the appropriate 
translation between the gravitational variables and quantities in the boundary theory. 
This dictionary allows us to write expressions for the entropy production and the 
change in other thermodynamic quantities induced by the quench. Section 7 provides 
a brief description of our numerical procedure. In section 8, we provide an independent 
analysis of the response in the slow quench limit, which later provides a check of our 
numerical results. We present and discuss various aspects of our numerical solutions 
in section 9. We conclude with a summary of the results and further comments in 
section 10. Finally, there are various appendices describing certain technical details. 
Appendix A presents explicit coefficients for the leading terms in the asymptotic of the 
expansion in section 3. Similarly, appendix B presents coefficients for the asymptotic 
expansions appearing in section 4. Finally, appendix C describes the variations of the 
renormalized bulk action constructed in section 5, which yield the expectation value of 
the stress tensor and the quenched operator. 
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2 Holographic model 



We will apply holographic techniques to study quantum quenches in a strongly coupled 
four-dimensional QFT. The quantum quench is implemented by adding a relevant 
operator with time-dependent coupling to the Lagrangian of the QFT, as follows [12]: 



In our calculations, the theory described by Lq is in fact a conformal field theory. The 
operator is relevant, meaning that it has conformal dimension A < 4. We will only 
consider A > 2 here, as natural in our holographic framework — see below. In our 
analysis, we start with the theory in a thermal state with A = and quench the system 
by switching on the coupling to some non-zero value. To further simplify our analysis, 
we will focus on quenches in the high temperature regime, where the temperature T 
provides the dominant scale in the problem. That is, we will only study quenches 
where A ^ stages. Note, however, that wc will allow the rate of change of 

coupling A to be arbitrarily large. In particular, we allow dt\{t) > T^~^. 

As the unperturbed QFT is a four-dimensional conformal field theory, the gravi- 
tational dual of the vacuum state is five-dimensional anti-dc Sitter (AdSs) spacetime. 
Since wc arc interested instead in a thermal state of the boundary CFT, the appropri- 
ate dual spacetime is an asymptotically AdSs planar black hole [13] — wc consider the 
boundary QFT in 7?^'^, hence the 'planar' geometry for the black hole horizon. Switch- 
ing on the coupling A is dual to switching on a massive scalar field in the gravitational 
theory. More precisely, we arc modifying the asymptotic boundary conditions for the 
bulk scalar field in a way that matches the profile \{t). Using holographic methods, we 
can easily determine the response of the QFT by examining the response of the scalar, 
which yields (Oa), as well as the response of the spacetime metric, which yields the 
energy density, pressure and entropy density of the boundary field theory. Since we are 
considering the high temperature regime, our calculations will be perturbative in the 
amplitude of the bulk scalar. That is, the scalar will only produce 'small' perturbations 
on the AdSs black hole background. 

The dimension A of the operator Oa is related to the AdS length scale L and the 
mass m of the bulk scalar field by [14] 



Notice that a relevant operator is dual to a scalar field with < 0. Of course, such a 



£0 ^ >Co + X{t) Oa. 



(2.1) 



A = 2 + V4 + L2m2. 



(2.2) 
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tachyonic mass is still consistent in five- dimensional AdS space as long as it respects the 
Breitenlohner-Freedman bound [15], i.e., m? > — In eq. (2.2), this imposes the 
constraint A > 2. The unitarity bound for a scalar operator in the four- dimensional 
CFT allows for A > 1, however, to study operators in the range 2 > A > 1, we must 
use the 'alternative quantization' of the dual bulk scalar set forward in [16]. However, 
we will not consider this possibility in the following and restrict our attention to A > 2. 

The dual gravitational theory is Einstein gravity coupled to a cosmological constant 
and a massive scalar field, i.e., 



J (fx^ (^i? + 12 - ^ {d^ f - im'^'^ . (2.3) 



_ 1 /■ , 1 ....9 1 



Since Newton's constant appears in an overall factor in front of the action, the scalar 
field $ is dimensionless. Further, we have also implicitly set the AdS curvature scale 
to one, i.e., L = 1, as can be inferred from the cosmological constant term. With this 
convention, it follows that irP' will also be a dimensionless number, which is implicitly 
given in units of 1/L^. As explained in [12], the scalar field might have further interac- 
tions, e.g., a <1>'^ potential, but any such higher order terms will not play a role in the 
following analysis describing the high temperature regime. 

As an aside, let us comment that it is natural to think of the unperturbed boundary 
theory as the A/" = 4 super- Yang-Mills (SYM) theory, in the limit of large N^. and large 
't Hooft coupling. In this case, our conventions are such that the five- dimensional 
Newton's constant is given by 

C 

However, we are slightly liberal in our analysis here in that we allow the conformal 
dimensions of Oa to take arbitrary values, rather than restricting ourselves to the 
spectrum of = 4 SYM. In this more general context, we can relate Newton's constant 
to the central charge of the boundary CFT with 

(2.5) 

where Ct is the central charge defining the leading singularity in the two-point corre- 
lator of the stress tensor — e.g., see [17]. 

Varying the action (2.3) with respect to the metric g and the scalar field one 
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obtains respectively Einstein's equations and the curved-space Klein-Gordon equation 
Q = E^, = R^, - ^ - g^, (^^ + 6 - J {d^f - , (2.6) 

= -^d^{,/^g''''d,^)-m^^. (2.7) 

We express our metric ansatz using infalling Eddington-Finkelstein (EF) coordinates 

ds^ = -A{v, y) dv^ + E'^{v, y) dx^ + 2 dvdy , (2.8) 

as was used in [9, 12, 18] in the context of holographic thermal systems. For the scalar 
in this background, we take $ = $(w, |/) (i.e., it is independent of the spatial directions 
x*). This choice allows us to describe homogeneous quenches where the coupling A 
is spatially constant but varies in time. The above is a convenient gauge (2.8) for 
numerically evolving the scalar field within a characteristic formulation. The resulting 
radial vector ^ is null and all points on a line with constant v (and x*) are causally 
connected. The resulting system of (partial differential) equations provide a nested 
system of (with both radial and time integrations) that can be evolved the spacetime 
radially from the boundary aX y = oo inwards and in forward in time. We will return 
to this discussion when we describe the numerics in section 7. 

With this metric ansatz (2.8), the Klein-Gordon equation (2.7) and Einstein's equa- 
tions (2.6) become [12] 

= 2E92/6) + 3 (9j^)$ + 3E9j/$ - m^E$ , (2.9) 
= i:d4t) + 2td^-2Y? + ^rn^^'^Y? , (2.10) 

= A + dlA-^td^ + i^dy^-^m''^\ (2.11) 

= E-isa^+^E($)2, (2.12) 

= 9^^E + ^E (9,$)^ (2.13) 
where we have defined for any function h{v,y), 

h = d^h+^Adyh. (2.14) 
More precisely, the above equations are obtained as: 

• Eq. (2.9) is equivalent to the Klein-Gordon equation (2.7) multiplied by E. 
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• Eq. (2.10) corresponds to the combination 



3 ^ 6 • 



yy 



0. 



(2.15) 



• Eq. (2.11) corresponds to the combination 



1 



{6Eii - SJ:''E^y - 4:AJ:''Eyy) = . 



(2.16) 



3S2 



Note that En denotes one of the diagonal components of Ejj^^, with /j, = v = i, 
i.e., there is no imphcit sum over i in this expression. 

• Eq. (2.12) corresponds to the combination 



• Eq. (2.13) corresponds to ^ Eyy — 0. 

Note that eqs. (2.12) and (2.13) are constraint equations, imphed by the previous three 
equations [12]. 

3 Solutions to the equations 
3.1 Static solutions 

As noted above, because we study quenches of the boundary QFT from an initial 
thermal state, wc consider the dual AdS spacetime initially containing a black hole. 
With $ = 0, the spacetime will have the static solution 



where the black hole horizon is located ai y — /j, and the asymptotic boundary of the 
spacetime is located at y = oo. This black hole solution gives the gravity descrip- 
tion of the original (conformal) boundary theory in thermal equilibrium. The QFT 
temperature is given by the temperature of the black hole, namely T — /i/tt.^ 

Now following [12], our analysis will be limited to considering a high temperature 
regime, where X{t) <S j-^-a noted above, this means that our calculations in the 

^Our conventions below will introduce a small correction to this result - see section 6. 



— — —ATlEyy — -:^A?YlEyy — . 



(2.17) 



^iv^y) = y^ — ^, 
y 

^v,y) = y, 



(3.1) 
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dual gravitational description are perturbative in the amplitude of the bulk scalar. In 
other words, we assume that the AdS spacetime contains a 'large' black hole and the 
scalar only makes 'small' perturbations on this background geometry. If we parameter- 
ize the amplitude of scalar field by the small parameter £, if follows from the Einstein 
equations (2.6) that the scalar only backreacts on the metric at order At the lowest 
order in the scalar and the metric can therefore be written as [12] 

A{v,y) = y^-t^ + ^^eA^{v,y) + o{e') , (3.2) 

y 



y2 

J:{v,y) = y + f,e'i:^{v,y) + o{f) , 



where factors of /i were introduced above to make both metric functions, Ap{v, y) and 
Ep(v,y), dimensionless. 

As a matter of convenience, we now change to the dimensionless coordinates p = 
///y, T = pLV, as well as x' = //x. For this choice of radial coordinate, the boundary 
hes at p = and the black hole horizon lies at p = 1. The scalar field and the metric 
coefficients are then written as 

$(r,p) = ^%{r,p) + o{e) , 

A{t,p) = p^[p-^-p^ + eAp{T,p) + o{^')) , (3.3) 
E(t,p) = /i(p-^+£%(T,p)+o(£^)) . 

In these coordinates, the metric then becomes 

ds'' = (-A(r, p) dr^ + E2(r, p) dx'^) - 2^ . (3.4) 

Note that the factor of pr'^ cancels with the p^ contained in the metric coefficients A 
and E^. The metric, and therefore the equations of motion will be independent of the 
black hole mass parameter p in these coordinates. 

If we consider the Klein-Gordon equation (2.9) to order I, the field $p decouples 
from the metric functions A^ and Ep and we are left with the linearized equation [12] 

+ 39,$p - (3 + p") 9,$p - 2pdrdp% + (p - p5) = 0. (3.5) 

The metric perturbations can then be determined from eqs. (2.10) and (2.11) at order 
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e [12]: 



= [-2(3-p^)+p2(i-p^)a2 + p(4a,-4a,-2pa,a,)]Sp 



m 



(3.6) 



24 



a.--(l-p^)(l + p9,) 



Ep + 2 [6 - 2pdp - p'^g 



777 

+ [2a.$p - (1 - p^) d,%\ + — $J . 



(3.7) 



Again, note that the mass parameter p does not appear in these equations (3.5)-(3.7). 

In the case of a static or equihbrium configuration, eq. (3.5) can be solved for the 
leading order scalar field 



%{P) 



4-A 4-A 4-A 



4 ' 4 

A\2 



-Cl 



r(^)r(|) 

fr(f) 



P^2Fl 



AAA 



(3.8) 



where 2F1 denotes a hypcrgeometric function. The constant Ci is arbitrary but the 
coefficient of the second term above is chosen to ensure regularity of the scalar at 
the horizon. Separately, both 2F1 (^^, p^) and 2F1 ( y, p^) have a 

logarithmic divergence near p = 1 but with the relative factor above, these logarithmic 
terms cancel in cq. (3.8). This static bulk solution will describe the system (to leading 
order in £) after it has equilibrated after the quench with a finite coupling A. Hence 
it will be useful to extract the relative magnitude of the normalizable and the non- 
normalizable modes of the bulk scalar in this new equilibrium configuration — see the 
next section. 



3.2 Time-dependent solutions 

In this subsection, we write down the asymptotic expansion for the leading order scalar 
$p(t, p) and metric functions, 74p(T, p) and Ep(T, p), in a time-dependent solution. 
Note that when A e Z or A e ^n+i {^-Q-j A = 2 or 3 as in [12]), logarithmic terms 
appear in these asymptotic expansions. However, generically these expansions do not 
contain any logarithmic terms and this is the case that we consider in the following. 
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The time-dependent solution $p(t, p) has an asymptotic expansion close to p = 
of the form: 



P'-^ (^(o)(r) + P0(o) + ^^^^^(0) + ^^'^'(0) + o (P^)) (3.9) 

Af^ /N 1 (2A-l)p2 (2A + l)p3. \ 

+P (^0(2A-4)(r) + p(/'(2A-4) + 4(A-l) '^^'^-^^ + i2(A-l) + ° )J ' 

where the coefficients 0(o) and 0(2A-4) are now functions of r. Here h = 9^-/1, for any r- 
dependent function h. In the following, we will choose some function for the coefficient 
of the non-normalizable mode, 0(o)(t), and then the normalizable coefficient 0(2A-4)(^) 
is determined by numerically integrating eq. (3.5). However, from the static solution 
(3.8), we have an analytic solution 

equihbrium : 0(2A-4) = A 2^. /^ "^(o) (^.10) 

for the late-time configuration describing the boundary theory after it has equilibrated 
with finite A. 

The solutions for the metric perturbations at order £^ take the form 

Mr,p) = 5^[a2,n(r)p"-^ + a2,„(r)/-^^+" + ^2,„(r)/^-^+"] , (3.11) 

n=4 

Sp(r,p) = ^[52,n(r)p"-^ + <72,n(r)/-^^+- + e2,„(r)p^^-^+"] , (3.12) 

n=5 

where (most of) the coefficients can be determined by solving eqs. (3.6) and (3.7) order 
by order in powers of p. However, the coefficient 02,4 enters these equations as a free 
parameter. Now taking the limit p — )■ 0, we simplify eq. (2.12) using results for the 
expansion coefficients from the other equations of motion to produce the following 
constraint: 

02,4 = ^ (a (2A - 5) (/)(2A-4)0(o) - (4 - A) (2A - 3) 0(o)0(2A-4)) , 

(3.13) 

and hence 

a2,4(r) = C - ^ (4 - A) (2A - 3) 0(0) (r) 0(2A-4)(r) 

+ ^ (A - 2) r dr' </)(2A-4)(t') 0(o)(t') , (3.14) 
J-00 
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where C is an integration constant. Following [12], we will choose C at a later stage 
so that the entropy production in the quench is proportional to 02,4 (t = 00). Note 
that since initially we have 0(o)(t = —00) = = 0(2A-4)(t = —00), it follows that 
a2,4(— 00) — C. Further if we set 0(o)(7" = 00) = 1, then 02,4 asymptotes to 

02,4(00) = a2,4(-oo) - ^ (4 - A) (2A - 3) 0(2A-4)(oo) 

2 1"°° 
+- (A - 2) / dr' 0(2A-4)(r') 0(o)(r') . (3.15) 

All the remaining coefficients appearing in eqs. (3.11) and (3.12) can be determined in 
terms of 0(o), 0(2A-4) and 02,4. Explicit expressions of some of the leading coefficients 
are given in appendix A. 



4 Fefferman-Graham coordinates 

We would hke to evaluate the entropy density, the expectation value of the stress-energy 
tensor and of the operator Ca in the boundary theory during a quench. Following the 
standard approach [19-21], we need to vary the on-shell gravitational action (2.3) with 
respect to the asymptotic boundary value of the appropriate fields — see section C. 
While EF coordinates are useful for evaluating the equations of motion, they are not as 
useful for determining the boundary one-point functions. The reason for the latter is 
that the "radial" direction dp is not orthogonal to the spacetime boundary located at 
p = 0, which is clear from the fact that the metric has off-diagonal r and p components. 
It will therefore be useful to transform to Fefferman-Graham (FG) coordinates [22] , in 
which the radial coordinate is orthogonal to the boundary of the spacetime. The FG 
coordinates have a spacehke radial coordinate r in contrast to the EF coordinates, with 
the null radial coordinate p. The FG coordinates are more appropriate for holographic 
renormalization, since we can choose a planar cut-off surface by simply fixing r to some 
small parameter e. 

In FG coordinates, the (asymptotically) AdS spacetime has the hne-element 

2 Gab{x,r)dx''dx^ dr^ 
= -2 + 7^' (4-1) 

a and b running from to 3. By equating this FG line-element (4.1) to the previous 
EF line-element (3.4) and writing the Eddington-Finkelstein coordinates r and p as 
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functions of the Fefferman- Graham coordinates t and r, we obtain a set of three equa- 
tions from which we can solve for r{t,r) and p{t,r), as well as the metric component 
Goo- The set of equations is 

= /x-MpVr' + (pr' + p't) , (4.2) 
-1 = {^^-'A {r'f + ^^p'r'^ , (4.3) 



Goo = ( -p-'Af' --fp] , (4.4) 



2 

where primes denote dr and dots denote df. We solve eqs. (4.2) and (4.3) by writing r 
and p as power series in r, with ^-dependent coefficients: 



^ n=l 

E ^(«) (^)^" + "-''"^ E (^)^" + E ^(«) (^)^^l ' (4-5) 

\n=5 n=0 n=l / 

p(t,r) = /ir + J]p(„)(i)r" + 

n=l 

f E (^)^" + ^'"'^ E (^)^" + ^'"^ E (^)^") ■ (4-6) 

\n=5 n=0 n=l / 

Upon solving for the above, we can also determine the metric Gab and scalar field $p 
in terms of similar asymptotic expansions in r 

( E C(n)a6(t)r" + r'-'^ Yl din)ab{t)r- + Yl Hn)ab{t)r A (,4.7) 

\n=4 n=0 n=0 / 

*p(^,r) = (/-'"^/(n)(t)r" + r^$^^(„)(t)r") . (4.8) 



n=0 n=0 



Explicit expressions of the leading coefficients are given in appendix B. For an asymp- 
totic solution of the nonlinear equations of motion in FG coordinates, see [23]. 

5 Holographic renormalization 

Given the metric and scalar field written in FG coordinates, we must evaluate the 
on-shell gravitational action (2.3). However, a naive evaluation yields a number of 
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divergences associated with integrating out to the asymptotic boundary at r = 0. 
Hence following the standard approach [19-21], we first regulate the calculation by 
introducing a cut-off surface r = e and then the divergences are eliminated by adding 
boundary countcrterms. Actually these counterterms are added in addition to the 
usual Gibbons-Hawking-Brown-York term 



Sghby = TTv / d^x^f^K 



(5.1) 



where 7a;) (e) is the induced metric on the cut-off surface and K is the trace of the 
extrinsic curvature of this surface. Recall that in our study, we choose the boundary 
geometry to be flat, i.e., 

g^^ = hmGab{t,r)=7jab, (5.2) 



and so the counterterm action turns out to be 



^ count 



4 - A 

d*xy/^\ -6 —^^ (5.3) 

where R (7) corresponds to the Ricci scalar constructed with ^ab- The (9$)^ and 
R (7) terms only cancel divergences which occur when A > 3 and so they should 
be discarded when A < 3. Although the term with R (7) vanishes to leading order 
when evaluated on a planar cut-off surface, it is required to cancel a divergence that 
arises in varying the metric to determine the stress tensor [24] . In particular, it cancels 
a divergent contribution to the pressure for A > 3 at order i'^. Also note that 
for the special cases A = 2, 3 and 4, there are also further logarithmic and finite 
counterterms, but we do not concern ourselves with these here. The interested reader 
can find a complete discussion of these cases in [12,25]. 

The holographic action Sreg — Sbuik + Sghby + S count can now be used to calculate 
the one-point correlators of the stress tensor and operator Oa- In order to calculate 
these expectation values, we need to vary Sj-eg with respect to the boundary metric and 
the scalar field, respectively. The details of these calculations are given in appendix C 
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and the final results are: 

SttG^S = _ ^2^4 0^^^^ ^ 1 ^2A - 3) (4 - A) 0(o)0(2A-4)) , (5.4) 

SnG^^^V = ^//^ - Qa2,4 - ^ (4A - 9) (4 - A) (/.(o)0(2A-4)) , (5.5) 
IGnCPiOA) = 2/x^£a;A(A-2)</.(2A-4). (5.6) 

Here S and V denote the energy density and pressure in the boundary theory, i.e., 
(T°°) — S and (T*-^) = 6^^ V. Further, ax is a proportionality constant relating the 
leading coefficient in the expansion (4.8) of the bulk scalar with the coupling in the 
boundary theory, i.e., £/(o) = ax A. We fix the precise value of this constant in section 
6.1 — see eq. (6.16). 

These one-point correlators must respect certain Ward identities [21]. In particular, 
one has the diffeomorphism Ward identity 

9^(r,,) = (e»A)9,A, (5.7) 

Of course, when the coupling A is constant, this expression reduces to the conservation 
of energy and momentum in the boundary theory. In the present case with a time- 
dependent coupling, the j — t component of eq. (5.7) yields 

dtS = -{OA) dtX. (5.8) 

Here the expression on the right-hand side describes the work done by varying the 
coupling in the boundary theory.^ Let us verify that eqs. (5.4) and (5.6) satisfy this 
constraint: First, comparing the expansions of the bulk scalar in eqs. (3.9) and (4.8) 
and recalling the relation £/(o) — axX from appendix C, we find to leading order 

0(0) =/^^"^aA^- (5.9) 
Then differentiating eq.(5.4), we find 

STTG^dtS = £V(-^a2,4-^(2A-3)(4-A)(0(o)0(2A-4)+0(O)0(2A-4))) 

= -£V(A-2) (/)(2A-4)0(o), (5.10) 

^Note that a minus sign appears here in accord with our conventions, which differ sUghtly from 
those in [12]. 
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where we simplified the expression by substituting for 02,4 from eq. (3.13). Now us- 
ing eqs. (5.6) and (5.9), we see that this expression precisely matches the expected 
Ward identity (5.8). Let us comment that this match should be no surprise since the 
constraint (2.12) (which was used to derive eq. (3.13)) reduces to precisely this Ward 
identity (5.8) on the asymptotic boundary r = [12]. 
We also have the conformal Ward identify 

T\ = {A-A) (Oa)A, (5.11) 

which follows from taking the trace of the stress-energy tensor with eqs. (5.4) and (5.5) 
and substituting eqs. (5.6) and (5.9). Here we do not find any anomalous terms (at 
quadratic order in i), since we are assuming that the operator Oa has a fractional 
conformal dimension. This result can be contrasted with the discussion in [12] which 
considered A = 2 and 3. 

6 Temperature and entropy density 

In this section we will calculate the temperature of the boundary theory before and 
after the quench, as well as the entropy produced during the quench. As described 
above, we are assuming that the quench takes the scalar field from a vanishing initial 
value with 0(o) = and 0(2a-4) = 0) to a final equilibrium solution where 0(o) = 1 
and 0(2A-4) = 0(2A-4)(oo)- In section 6.3, we will consider 'reverse' quenches which 
instead take the system from 0(o) = 1 to 0. In our perturbative calculations for high 
temperature quenches, we find that if the profile for the 'reverse' quench is given by 
0(0) (r) — 1 — 0(0) (r), where 0(o)(t) describes some 'forward' quench, then we find 
that 0(2A-4)('?') = 0(2A-4)(oo) — 0(2A-4) (t) , where (j)(2A-4){T) is the response for the 
corresponding 'forward' quench. Similarly, we will find that the entropy production 
is the same in the forward and reverse quenches. Further, in the case of an adiabatic 
quench, no entropy is created and the process is reversible. 

As discussed in section 3.1, the initial configuration before the quench is the well- 
known planar AdS black hole described by eq. (3.1). The calculation of the corre- 
sponding temperature is a straightforward exercise with the result T — /i/tt. However, 
recall that in eq. (3.14) we established a convention where a2,4(— 00) = C. That is, our 
metric perturbation is nonvanishing even at r = —00. The effect of this convention 
is to shift the black hole mass parameter, i.e., /i — >■ /i^ where = 1 — £^02,4 (—00). 
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Hence, to quadratic order in the expansion in the initial temperature becomes 



6.1 Final temperature 

Next we wish to determine the final equilibrium temperature of the system after the 
quench has taken place. This calculation is more subtle as with our perturbative calcu- 
lations, since we will not have the full metric describing the final black hole geometry. 
Instead then, we turn to the thermodynamics of the boundary theory to determine the 
final temperature. That is, we will compare the energy density and pressure in QFT 
variables (already in terms of the final temperature Tf and the coupling A) to the energy 
density and pressure calculated holographically in terms of gravitational variables. In 
doing so, we are able to derive meaningful relations between the field theory coupling 
and temperature and the bulk parameters // and I. Of course, by assuming a form for 
E and V-i our final temperature and entropy production will necessarily depend on the 
conventions used to define our coupling. This cannot be helped, because we do not 
know the Lagrangian for the boundary theory when the quench is by an operator of 
arbitrary dimension A. This can be contrasted with the discussion in [12] for the cases 
of A = 2, 3, where the exact equilibrium expressions for E and V are known from [26]. 
Nonetheless, we will find physically meaningful interpretations for our results. 

To begin, we make the following ansatz for the energy density and pressure in the 
final equilibrium of the boundary theory. 



where \f — A(t = oo) denotes the final value of the coupling. To leading order our 
ansatz reduces to the expressions expected for a conformal theory and is in accord 
with our analysis, the perturbation of these conformal terms is quadratic in the cou- 
pling. Further, we have expressed the perturbations in terms of the dimensionless ratio 



normalization for the coupling. We can compare these expressions with those given 




(6.1) 





(6.3) 



Xf/Tf-^. Setting the pre-factor for this term in the pressure (6.3) really defines our 
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in [12,26]. For example, we find for A = 3, 



2r(|) 



4 



Aj^^mJ, (6.4) 

where ruf was tfie fermion mass in tlie boundary theory. Using this expression, we can 
confirm the results derived below for the equilibrium values 
of the observables agree with those given in [12,26]. 

Now we need to determine the constant of proportionality a/ in eq. (6.3). To 
proceed, we only assume that the boundary theory obeys standard thermodynamics, 
following [27]. First, we write the free energy density as 

F^E-TS, (6.5) 

where S is the entropy density. In the absence of any chemical potentials, F = —V. 
Therefore combining these expressions with eqs. (6.2) and (6.3), the final entropy den- 
sity is given by 

A 




3^/ |4-(3«/ + 1)|;m^) I ■ (6-6) 

We use the first law of thermodynamics (with fixed volume) to write 

dSf _ dS 



dTf ' dTf 



The left-hand side of eq. (6.7) is 



dS 



whereas the right-hand side is 




T,^ = ^rH4-i(3a, + l) (2A-5) 




By comparing these two expressions, we solve for as 

2A - 5 , , 

«/ = ■ (6.8) 

Note that it may seem that the quench has no effect on the energy density for A = | 
(when af = 0), but even in this case, the initial and final temperatures will differ by a 
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term of order Xj. Hence, there will still be a change in S in this case, contained in the 
term in eq. (6.2). 

Next, we compare these results for the boundary theory with the corresponding 
expression in the gravitational dual. In particular, we would like to find £ in terms 
of the temperature Tf and the coupling A/. However, first we fix the normalization 
factor A appearing in eqs. (6.2) and (6.3). This factor would be the unchanged in the 
initial equilibrium of the conformal boundary theory, i.e., at t = — oo, we would have 
£i — AT^. Comparing the latter expression with eq. (5.4) then yields 

ATt^ / (1 - f a2,4(-oo)) . (6.9) 

Given the expression for the initial temperature in eq. (6.1), we see that 

Next, we take the trace of the stress tensor in both the field theory and the gravi- 
tational dual: 

2 



(Tqpt)\ = -^^T;(4-A)(^) , (6.11) 



(TgrT a = (4 - A) (A - 2) 0(o)0(2A-4) 



.(5) 

where 0(2A-4)(oo) is given by eq. (3.10) with 0(o) = 1, i.e., 



(4-A)(A-2)0(2A-4)(oo), (6.12) 



= (6-13) 

Equating the two expressions above and using eq. (6.10), we find 

to leading order in \f/T^~^. Note that here we have also used eq. (6.1) to substitute 
/x^ = vr^Tj + o(Aj) since the initial and final temperatures will only differ by o(Ay.) in 
our perturbative calculations. Further, the above expression takes account of the fact 
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that 0(2A-4)(oo) is always negative in the range of interest, i.e., 2 < A < 4 — see 
eq. (6.13) above. RccalUng that we set 0(o)(c)o) = 1, we note that imphcitly the right- 
hand side of eq. (6.14) is actually ^^0(o) ^^'■^ ^'^ ^^^^ equation fixes the normalization 
between the leading coefficient in the asymptotic expansion of the bulk scalar and the 
boundary coupling, i.e., 

^ +o{X'). (6.15) 



V(A-2)|0(2A-4)(oo)| T^-^ 

Alternatively in appendix C, we introduced the proportionality constant ax in i'/(o) = 
ax A. So comparing the expansions of the bulk scalar in eqs. (3.9) and (4.8) using eq. 
(B.37) , we now have 

^ /(A otu ( =T\ + ^^^'^ ' ^^-^^^ 

^(A-2) |(/>(2A-4)(00)| 

where as above, we used /i^ — n^T^ + o(A^). 

6.2 Entropy production during the quench 

Here we extend the previous analysis to determine the entropy production during the 
quench. First using the expression for the free energy density (6.5), as well as F = —V, 
we find 

.q, T. P.. -^ -p. 

(6.17) 



Si Tf £i + Vi 

Initially the boundary theory is conformal and the vanishing trace of the stress tensor 
requires Ei — 3Vi. Now the latter can be used to re-express eq. (6.17) as 



m + . (6.18) 

First, we determine the ratio of the temperatures by equating the final energy 
densities given in terms of the gravitational variables (5.4) and of the boundary theory 
(6.2). The initial temperature is introduced here by substituting for /i using eq. (6.1), 
which then yields 

^ = 1 + (^02,4(00) - a2,4(-oo) + ^ (2A2 - 8A + 9) 0(2A-4)(oo)^ . (6.19) 
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Now using the expressions for the energy density and pressure in eqs. (5.4) and (5.5) 
at the initial and final times, we find: 



Si 
Vi 



1 - (^a2,4(oo) - a2,4(-oo) + ^ (2A - 3) (4 - A) 0(2A-4)(oo)^ , (6.20) 
1 - (^02,4(00) - a2,4(-oo) - ^ (4A - 9) (4 - A) 0(2A-4)(oo)^ . (6.21) 



Combining these results in eq. (6.18) then yields 

3£2 



Si 
Si 



1 - 



^02,4(00) - a2,4(-oo) - ^ (A - 3) (A - 1) 0(2A-4)(oo)^ . (6. 



22) 



Now recall from eq. (3.14) that a2,4(— 00) = C, where the latter is an arbitrary inte- 
gration constant. Hence following [12], we choose this constant to simplify the above 
ratio of entropies, i.e., 



a2,4(-00) = -- (A - 3) (A - 1) 0(2A-4)(OO) . 



(6.23) 



Hence, after substituting for and a2,4(— c>o) from eqs. (6.14) and (6.23), respec- 
tively, the ratio of the final and initial entropies (6.22) becomes 



5/ _ _^ 3 02,4(00) / 

Si ~ ^4(A-2)0(2A-4)(oo) It; 



/ 



4-A 



(6.24) 



Further substituting for i'^ and a2,4(— 00) in eqs. (6.19)-(6.21), we find the change in 
temperature, energy density and pressure are given by 

2 



AT 

Ti 


"A -2 
[ 6 + 


A£ 


" 1 

_ 3 ^ (A ■ 


AV 


"2A-7 


Vi 


3 



02,4(00) 



6 4 (A-2)0(2A-4)(oo) 



A 



/ 



rp4-A 



02,4(00 



+ 



(2A-4)(oo) 
02,4(00) 



A/ 



t; 



4-A 



(A -2) 0(2A-4)(OO) 



A/ 



rpi-A 



(6.25) 
(6.26) 
(6.27) 



where our notation is e.g., AV = Vf — Vi. 

The second law of thermodynamics demands that the ratio Sf/Si must always be 
greater than one. Hence requiring 02,4(00) < becomes a test of our numerical solutions 
and we successfully confirm that this inequality is satisfied in the obtained numerical 
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solutions. Since </>(2A-4)(c>o) is always negative in our analysis (and we restrict our 
attention to A > 2), eqs. (6.25) and (6.26) indicate that the changes in the temperature 
and the energy density are always positive. However, from eq. (6.27), the change in 
pressure is only guaranteed to be positive for A > 7/2. Otherwise, the pressure can 
either increase or decrease depending on the precise value of A and the magnitude of 
02,4(00). A more detailed discussion is given in section 9.5, where we consider the effect 
of the numerically determined values of 02,4(00) on the shifts of these quantities. 

Another check of the present analysis comes from considering the adiabatic limit. 
As we discuss in section 8 in this case, the system remains in a quasi-static equilibrium 
with 0(2A-4)(^) = 0(2A-4)(oo) 0(o)(t). Substituting this expression into eq. (3.15), as 
well as using the integration constant chosen in eq. (6.23), it is straightforward to show 
that 02,4(00) vanishes. Hence as expected for an adiabatic transition, no entropy is 
produced, as discussed in [12]. 

As a final consistency check, we consider the speed of sound in the thermal plasma, 
which is given by 

" " [dTfJ I [dTfJ 

- ^^(4-A)(A-2)(^) . (6.28) 

Note the second term is negative for all A in the range 2 < A < 4. Hence we find 
c1 < 1/3, as required by [28]. While = 1/3 for A = 2 and 4, our analysis only 
applies for the conformal dimension strictly limited within the range 2 < A < 4. 

6.3 Reverse quenches 

Up until now we have assumed that the quenches begin with the boundary theory being 
conformal, i.e., X — and then end with some finite A. In the gravitational description 
then, they involve some profile 0(o)('?') which begins with 0(o) = at r = —00 and ends 
with 0(0) = 1 at T = 00. In this section, we consider 'reverse' quenches in which the 
coupling is initially finite and is brought down to zero. In particular, we can readily 
repeat the analysis for reverse quenches where the non-normalizable coefficient of the 
bulks scalar is chosen to be 

0(o)(t) = 1-0(o)(t). (6.29) 
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Because our analysis is limited to the perturbative high temperature regime, the equa- 
tion of motion (3.5) for the bulk scalar is linear and hence we can add any two solutions 
to produce a third solution. In particular then, adding the scalar field solutions for the 
forward and reverse quench must yield the equilibrium solution with 0(o)(t) — 1. Alter- 
natively, the reverse quench produced from eq. (6.29) is simply the equilibrium solution 
(3.8) (with Ci = 1) minus the time-dependent solution describing the forward quench. 
In the equilibrium case, the normalizable coefficient in the bulk scalar is 0(2A-4)(oo) 
and so the corresponding coefficient in the reverse quench must be 

0(2A-4)(t) = 0(2A-4)(OO) - <;i(2A-4) (t) , (6.30) 

where 0(2A-4)('?") denotes the response produced in the original (forward) quench. 

In the reverse quench, the metric coefficient 02,4 still satisfies eq. (3.14) and so we 
have the solution 

a2,4(T) = C-^(4-A)(2A-3)0(o)(t)0(2A-4)(t) 

+ ^ (A - 2) ciTV(2A-4)(T') ^(0)(t') . (6.31) 

where C is a new integration constant. Note that in this case, the second term van- 
ishes for r — 7> 00 but as r — > —00, it is proportional to 0(2A-4)(— 00) = 0(2A-4)(oo). 
Repeating the analysis of the previous section for our reverse quench and demanding 
that the entropy production is now proportional to 02,4, we find that the integration 
constant must be chosen as 

C = ^(A-2)0(2A-4)(-oo). (6.32) 
With this choice then, we have 

bi 4(A-2)0(2A-4)(-oo) j 

Further, when we compare the expression for 02,4(00) with eq. (3.15) for 02,4(00), it is 
straightforward to show that these two constants are equal, i.e., 

02,4(00) = 02,4(00) , (6.34) 

just as was found in [12]. Hence comparing to eqs. (6.24) and (6.33) and noting that 
0(2A-4)(— 00) = 0(2A-4)(oo), we scc that the entropy production is identical in the 
forward and reverse quenches. 
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We can also find the changes in the temperature, the energy density and pressure 
as before: 



AT 
A£ 
AV 



A - 2 1 

+ - 



02,4(00) 



6 



1 
3 



4 (A-2)0(2A-4)(OO 

02,4(00) 



(A -2) 0(2A-4)(OO) 



2A-7 



+ 



02,4(00) 



(A -2) 0(2A-4)(OO) 




(6.35) 
(6.36) 
(6.37) 



where again e.g., AV — Vf — Vi- Comparing these expressions for the reverse quench 
with the corresponding results in eqs. (6.25)-(6.27) for the forward case, we see that 
02,4(00) and 02,4(00) have the same coefficients while the constant terms are equal but 
with opposite sign. Hence for the reverse quenches, whether any of these physical 
quantities increases or decreases depends on both the magnitude of 02,4(00) and the 
value of A — see section 9.5 for a further discussion. We also note that in the case of 
an adiabatic quench (see section 8), we find 02,4(00) ~ and hence 02,4(00) = 0. That 
is, for adiabatic transitions, no entropy is produced. However, we would also find that 
the changes in the temperature, energy density and pressure are exactly opposite in 
the forward and reverse cases. 

Since the forward and reverse quenches are simply related in our perturbative high 
temperature analysis, we will continue to focus on the forward quenches in the rest of 
the paper. 



7 Numerical procedure 

We now briefiy describe the numerical simulations which we used to understand the 
quenches. Essentially we implemented the same approach as in [12], using numerical 
techniques developed in [29,30]. For details on our numerical implementation the in- 
terested reader can see the appendix in [12]. The primary purpose of our simulations 
was to find the response 0(2A-4)(t) for a given source 0(o)(t), as described earlier, it is 
convenient to work with infalling characteristics and in particular we adopt Eddington- 
Finkelstein coordinates as in eq. (2.8) or (3.4). These coordinates are regular at the 
horizon allowing us to excise the black hole from the computational domain (by inte- 
grating some distance inwards and stopping the integration as this region is causally 
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disconnected from the outside). 

We choose the source term in the asymptotic expansion (3.9) of the bulk scalar to 

be 

0(o)(r) = ^ + ^tanh(^). (7.1) 

This describes a family of quenches where, as desired, the source starts at zero in the 
asymptotic past and ends at one in the asymptotic future. Here a controls the rate at 
which the quench takes place. In terms of the dimensionful time, we have t — /iv /it 
and hence the timescale on which the transition from zero to finite coupling is made is 

At^-^—. (7.2) 

Hence for a ^ 1, the quenches are 'slow' which means that the transition occurs on a 
timescale that is much longer than the thermal relaxation timescale. Alternatively for 
q; ^ 1, the quenches are 'fast', meaning the transition timescale is much shorter than 
the thermal timescale. The limit a — > would correspond to an 'instantaneous' quench, 
i.e., 0(0) becomes a step-function. The limit a ^ oo corresponds to an adiabatic 
transition — see section 8. We note that with cq. (7.1), the source profile has a 
continuous derivative for all time. More general quenches which arc not infinitely 
differentiable with respect to time will be studied in a later paper [31]. 

Again the goal of our simulations is to find the response 0(2A-4)(t) for a given 
source 0(o)(t). Practically, it is more convenient to solve the linearized scalar equation 
(3.5) in terms of 0(r, p), which is defined by 

%{t:P) = P'-^ ^0(0) + P0(O) + ^^^^0(0) + ■■■)+ P'' 0(t, P) . (7.3) 

Here the term in brackets contains the leading terms in the asymptotic expansion (3.9) 
of $p. In particular, we include any terms which are leading compared to p^. Formally 
then the leading solution takes the form = P^~'^0(2A-4) + • • • and so the leading 
behaviour in contains the desired response function. However, as a practical matter, 
we obtained noticeably better accuracy for 0(2A-4) by fitting (at each timestep) with 
an expansion of the form 

0(t, p) = eop'-^-^ + eip^-^--^ + • • • + p^-^<P(2A-4) + Sy-^-^^'^ + o , (7.4) 

where the next leading order after p^ The coefficients are, of course, 

all very small since the terms at these orders are already included in eq. (7.3). The 
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factor 6n is typically not small, but including this term nonetheless gives better results 
when fitting 0(2A-4)- In the cases that we choose k > 4 — A, we only include terms in 
eq. (7.4) with nonnegative powers of p. 

We introduce k for convenience in the fit and choose it so that (f) still vanishes at 
the asymptotic boundary, i.e., p — 0. We further choose k so that the relative power of 
p between the two contributions in eq. (7.3) is an integer, i.e., A + n is an integer. This 
ensures that upon substituting into the equation of motion (3.5), after simplification, 
only integer powers of p appear in the coefficients of and its derivatives, as well as 
in the source terms. Our choices of k, for each of the conformal dimensions considered 
in our simulations, are shown in table 1. 

Table 1: Choices made for k while simulating 4^{t,p) for various A. 



A 


7/3 8/3 10/3 11/3 4 




5/3 1/3 5/3 1/3 2 



Once the numerical solution has been obtained at each timestep, we fit (f) with a 
series in p with rational exponents, as described above, to determine the coefficient 
0(2A-4)- Repeating this process for each timestep then generates the full profile for 

0(2A-4)(t). 



8 Slow quenches and the adiabatic limit 

In this section, we consider the case where the transition between the initial and final 
theories is made arbitrarily slow. In fact, we find an analytic solution of 02,4 for such 
slow quenches below. The results derived from this approach provide an independent 
check for our numerical solutions — see the discussion in section 10. 

Let us first consider the adiabatic limit: Given our choice for the integration con- 
stant a2,4(— 00) in eq. (6.23), the expression (3.15) for 02,4(00) after equilibration be- 
comes 

1 2 ■ 

a2,4(oo) = - 3 (A - 2) 0(2A-4)(oo) + - (A - 2) j drV(2A-4) (r') 0(o)(t') . (8.1) 

As noted in [12], in the adiabatic limit, the system remains in a quasi-static equilibrium 
throughout the transition. Hence the ratio of the normalizable and non-normalizable 
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coefficients in tlie bulk scalar are precisely as given in eq. (3.10) at every stage of the 
transition, i.e., 

adiabatic: 0(2A-4)(t) = 0(2A-4)(oo) 0(o)(r) . (8.2) 
In this case, the integral in eq. (8.1) becomes 

/ C^'r>(2A-4)(T')0(O)(T') = ^0(2A-4)(OO) / dr' 9^/ (0(o) (t')) 
J —oo J —oo 

= ^0{2A-4)(OO). (8.3) 

Given this result, we easily see that 02,4(0x3) (and hence the entropy production) van- 
ishes in the adiabatic limit. 

For arbitrarily slow quenches, the response should approach the adiabatic profile 
(8.2) as a limit. Given the profile in eq. (7.1), this slow limit is achieved by taking a 
large. Then following [12], we can write the bulk scalar $p in an expansion in inverse 
powers of a, i.e., 

00 

% = 0(0) (r) $e(p) + E 0(")(r) Ri"\p) , (8.4) 

n=l 

where 0(o) (r) is the profile in eq. (7.1) and $e(p) is the equilibrium solution given in 
eq. (3.8) (with ci — 1). Hence the first term above is the desired solution describing an 
adiabatic transition, i.e., this term yields the response in eq. (8.2). We have assumed 
that each of the subsequent terms in the series are separable and this form is easily 
confirmed, e.g., see below. 

Let us solve for the first correction in eq. (8.4), (j)^^\T) Ri^\p) , since this term 
produces the leading contribution to 02,4 in the slow quench limit. Substituting our 
ansatz (8.4) into the decoupled Klein-Gordon equation (3.5), we find at order a"^ 

^i^Hr) (1 - dl + p + 3) d, + m'] i?W(p) 

= a 0(0) p(3-2p9p)$e(p). (8.5) 

Note that a 0(0) is order a*^ in our expansion because the derivative of 0(o) brings out 
a factor of 1/a. Now solving eq. (8.5), requires that the time-dependent function takes 
the form 0^^^ (r) = a 0(o) while the radial profile satisfies the inhomogeneous equation: 

(1 - P') dl + p{p' + 3)d, + A{A- 4)] (p) = p (3 - 2p d,) $e(p) . (8.6) 



27 



The solution to the above equation can be written in series form as 

oo oo 

R'^^\p) = p'-^ E ^i-)P^ + P^Y. hn)P'' ■ (8.7) 

n=0 n=0 

Here we have two independent integration constants, a(o) and 6(o). The first coefficient 
a(o) is set to zero, since the order contribution to $p defines the source and we do 
not want any new sources beyond 0(o) to appear at higher orders in the 1/a expansion 
in eq. (8.4). To determine 6(o), we demand that Ri^\p) is regular on the event horizon 
at p = 1. To analyze the profile close to the horizon, we make the change of coordinates 
z — 1 — p and solve for Ri^\z). Near z — 0, we may write Ri^\z) as the series 

oo 

i^W(^) = Ec(„)^^ (8.8) 

n=0 

which includes only the solution that is well-behaved a.t z = 0. In this series, we now 
have the undetermined coefficient C(o). 

To proceed further, we resorted to solving eq. (8.6) numerically. In particular, we 
produced independent solutions integrating in from z = 1 (or p = 0) and integrating 
out from z = (or p = 1). Then by matching the two solutions for Rs^^ at an inter- 
mediate point between the asymptotic boundary and horizon, we solved for both 6(o) 
and C(o) simultaneously. This shooting method was used to solve for the undetermined 
coefficients in the cases A = |, |, y and We should emphasize that although this 
approach is again numerical, it is an independent approach very different in spirit from 
the numerical approach described previously.^ 

Comparing the present solution with the expansion in eq. (3.9), we see that the 
corrected response takes the form 

0(2A-4) (t) = 0(2A-4) (oo) 0(o) (t) - 6(0) 0(0) (t) + o(1/q;^) . (8.9) 

Again the second term above is of order 1/a because the derivative acting on 0(o)('7') 

brings out this factor. Now upon substituting this expression into eq. (8.1), the adia- 

batic response yields zero and so we are left with 

2 f°° 
02,4(00) = --(A-2)b^o)J_ c^TV(2)(r') + 0(1/^2) 

" -|^(^-2)^(o)+o(l/«'), (8.10) 

^It is straightforward to write a formal Green's function solution for cq. (8.6). Given this form of 
the solution, the coefficient 6(0-) can be determined by numerically evaluating a specific integral. The 
results produced this way agree well with those given in eq. (8.11). 
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where the final result is produced by inserting eq. (7.1) for 0(o)- Hence given the 
numerical solution for 6(o) for each value of A listed above, we can determine the 
leading contribution to 02,4 for slow quenches. Our results are as follows: 

A=-: 024(00) = -0.01958 -, 

3 ' a 

8 1 
A=-: 024(00) = -0.05205 -, 

3 ' a 

10 . . 1 

A = — : a2,4 (00) = -0.09838 - , 

o Qi 

A = ^: 024(00) = -0.1083 -. (8.11) 

3 ' a 

These results compare well to our numerical results, as discussed below in section 10. 



9 Results 

We now turn to the results of our numerical simulations. We determined the response 
functions for various values of a and with several different masses of the bulk scalar. 
Explicitly, the scalar masses for which we studied the quenches are: — ~Y ("^i^^ 
A = |), = -f (with A = |), = -f (with A = ^) and = The masses 
were chosen so that the conformal dimension lies in the desired range 2 < A < 4 and 
so that 4 — 2A is not an integer. The latter ensures that the asymptotic expansion 
(3.9) for the scalar does not contain any logarithmic terms. We further comment on 
A = 4 case, which together with results of [12] would form a complete picture of 
perturbativc quenches in strongly couple gauge theories induced by the coupling of a 
relevant operator. 

In subsection 9.1, we extract the response function d>(2A-A) for fast quenches {i.e., 
a < 1) from our numerical data for each A. From that we sec that 02,4 and the 
maximum displacement of 0(2A-4) follows a scaling behaviour for fast quenches. In 
subsection 9.2, we see that the scaling of these quantities follow a universal behaviour, 
determined by A and a only. In subsection 9.3, we show the numerical results for slow 
quenches {i.e., a > 1). We show that for very large a, 0(2A-4) approaches the adiabatic 
response (8.2). We also find that 02,4 scales as a'^, as expected from the analysis in 
section 8. Finally, in subsection 9.4, we study the excitation times of the response, 
by considering when it first deviates from the adiabatic response by more than 5% 
(an arbitrary threshold). We find that for fast quenches, the scaled excitation time 
Tex/c( scales logarithmically with the quenching parameter a. We also see a universal 
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behaviour in the slopes — y ^log a j quench limit for fractional A. In the 

same subsection we study the relaxation times for the quenches, defined by the last 
time that 0(2A-4) falls below a 5% deviation from its final equilibrium value. 

9.1 Response for fast quenches 

Here we hst the results for intermediate to fast quenches {i.e., a < 1). We first find the 
response 0(2a-4) and then determine the coefficient 02,4(00) using equation (3.15), as 
well as the integration constant a2,4(— 00) in eq. (6.23). By doing multiple simulations of 
the time-dependent behaviour of 0(2A-4) with increasingly finer discretizations, we can 
extrapolate for the continuum value of 02,4(00). Knowing 02,4(00), we can determine 
the response of various physical quantities in the boundary theory to the quench, using 
eqs. (6.24)-(6.27). Note that these results would also apply to the reverse quenches, 
as discussed in section 6.3. We also study the behaviour of 0(2A-4) and 02,4(00) as 
functions of a. In the next subsection, we compare this a dependence for different 
values of A and identify the universal behaviour for operators of different dimension. 

In figure 1, we plotted 0(2A-4) against time for various values of a. The time is 
rescaled by a factor of l/a so that the different plots equilibrate in approximately the 
same distance on the horizontal axis. In each of the plots, we have also rescaled the 
vertical axis by a^^~^. Hence the maximum displacement of 0(2A-4) is actually growing 
with decreasing a. With the rescaled axes, the peaks corresponding to smaller a lie 
close together, which seems to indicate that there is a universal scaling depending on 
the operator dimension. Further with these scalings, the response seems to converge to 
a particular limit with decreasing a. For fast enough quenches, we expect these plots 
to coincide exactly. 

In figure 2, we show a different visualization of the same fast quenches. Here, 
q;^^~'*0(2A-4) is plotted as a function of the source 0(o). It is perhaps easier to see 
the convergence of the faster quenches to a limiting curve. Note that because of the 
growth of 0(2A-4) for fast quenches, the expression (3.15) for 02,4(00) is dominated by 
the integral. Further the latter can be re-expressed as o?0(o)0(2A-4) and so for these 
fast quenches, 02,4(00) is essentially given by the area under the curves in figure 2. 

In figure 3,^ we see that log (—02,4(00)) plotted against logo; tends to a straight 
line for small values of a. This indicates that 02,4(00) scales as some power law of the 
quenching parameter a for very fast quenches. A fit of this linear behaviour suggests 

^Note that we plot the logarithm of —02,4(00) since 02,4(00) is always negative. 
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Figure 1: (Colour online) Plots of the response coefficient 0(2A-4) for quenches of 
different speeds, in the fast quench regime. Time is rescaled by a factor of l/a and 
the value of 0(2A-4) is rescaled by Clockwise from the top left, the plots are 

for A = 7/3, 8/3, 11/3 and 10/3. In each case, the response is presented for a = 1 
(dashed), 1/2 (brown), 1/4 (blue), 1/8 (purple), 1/16 (green), 1/32 (orange) and 1/64 
(red) (as well as a = 1/128 (yellow) for A = 7/3). 

the slope matches 4 — 2A in each case. The lines shown in the plot are the linear 
fits through the points corresponding to the three fastest quenches, thus showing the 
asymptotic behaviour of 02,4(00) for fast quenches. Although it is not shown, for 
small values of a, the logarithm of max{|0(2A-4)|} plotted against logo; also tends to 
a straight line with the same slope as in the plot of log(— 02,4(00)). 

As discussed in [12,31], the scaling of the response 0(2 a-4) for fast quenches is more 
subtle when A = 2, 3 and 4. Specifically, it is rather the 'subtracted' response, defined 
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Figure 2: (Colour online) Plots of 0(2A-4) for quenches of different speeds, in the fast 
quench regime. 0(2A-4) is rescaled by a^^"^ and is plotted against the actual value 
of the source 0(o). Clockwise from the top left, the plots are for A = 7/3, 8/3, 11/3 
and 10/3. In each case, the response is presented for various values of a, which are 
indicated using the same colour scheme as in figure 1. 
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(2A-4) 



^>2A-4 + < 
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Id 
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for A = 2 , 
for A = 3 , 
for A = 4 , 



which scales faithfully in the limit of fast quenches, i.e., 



{9.1] 



lim 



a-s>0 



(2A-4) 



constant . 



(9.2) 



Notice, however, that the additional In a terms in (/)2A-4 above do not contribute to 
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Iog(-a2,4(oo)) log(-a2,4(oo)) 




Figure 3: log-log plots for —02,4(00) versus a for various A, in the fast quench regime. 
The straight lines shown are least-squares fits through the three leftmost data-points 
in each case. The fact that the plots tend to straight lines for negative values of log a 
means that 02.4(00) scales as a power law for small a. Clockwise from the top left, the 
plots are for A = 7/3, 8/3, 11/3 and 10/3. 



02,4(00) for A = 3 or 4. For these cases, we have 



A = 4 



/oo 
^'^ <?^(o)(t) 0(0) (t) = 
■00 

/oo 
^'^ 0(O)('r)0(O)(T) = 
-00 



1 • 



0(0) (t)- 



0(0) ('t) 0(0) (t) - 7^0(0) (t) 



(9.3) 



= 



which, as we have indicated, both vanish for a generic source (/)(o) as long as the profile 
becomes constant as r — )■ ±00. Hence for A = 3 and 4, we still have 02,4(00) ~ l/a^^~^ 
for fast quenches as a — )■ 0. On the other hand, as shown in [12], the logarithmic term 
in eq. (9.1) gives the dominant contribution in eq. (3.15) for fast quenches with A = 2 
and we have instead 



02,4(00) = ^ \Vi0L + 0{\) , 



as q; — > . 



(9.4) 
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(ilog(-a2,4(oo)) 
dloga 
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Figure 4: (Colour online) Plot of the asymptotic scaling of —02,4(00) as a power of a, in 
the fast quench hmit. The line shown is the predicted theoretical trend = 
2A -4. Points shown are for A = 2, 7/3, 8/3, 3, 10/3, 11/3 and 4. The datapoints for 
A = 2, 3 are taken from [12]. The datapoint for A = 4 is taken from [31]. 

9.2 Universal behaviour for fast quenches 

For fast quenches, our numerics above suggested that —02,4(00) grows as 0;"^^^"^^ as 
a becomes arbitrarily small. This behaviour is confirmed in figure 4. This plot shows 
~''^°1roga'^°°" v^^^^s A for A = 2, A = 7/3, A = 8/3, A = 3, A = 10/3, A = 11/3 and 
A = 4. The data-points are the slopes of the straight line fits for the plots shown in 
figure 3. The data-points for A = 2 and A = 3 are taken from [12] while the data-point 
for A = 4 is taken from [31]. We set the log derivative to zero for A = 2 but, as noted 
above in eq. (9.4), 02,4(00) actually scales as a logarithm of a [12]. 
As well as the individual data-points, we have plotted the line 

-^M^^mM = 2A-4. (9.5) 
ologo; 

In fact, a least-squares fit through the data-points yields 

-^^^i=^MMl = i.99A_3.96, (9.6) 
o log a 

which matches the expected trend within our numerical accuracy. Hence figure 4 
confirms the universal scaling 

1^2,4(00)1 oc (9.7) 
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Figure 5: (Colour online) Plots of 0(2A-4) for quenches of different speeds in the slow 
quench regime. Time is rescaled by a factor of and the value of 0(2A-4) is rescaled 
by Plots for larger a follow an inverted tanh-profile more closely, which is 

a negative constant times the source. Clockwise from the top left, the plots are for 
A = 7/3, 8/3, 11/3 and 10/3. In each case, the response is presented for a = 1 
(dashed), 2 (brown), 4 (blue), 8 (purple), 16 (green), 32 (orange) and 64 (red). 

for fast quenches by an operator with 2 < A < 4. 

As noted above, the maximum displacement of </)(2A-4) seems to exhibit the same 
scaling behaviour as above when A is fractional. However, as also commented above for 
A = 2, both 02,4(0x3) and 0(2A-4) grow as — logo; for small a [12]. Similarly, for small 
a with A = 3 and 4, 02,4(0x3) has the above scaling but the maximum displacement of 
0(2A-4) exhibits an additional log a growth on top of this simple scaling — see further 
discussion around eq. (9.1) and in [31]. 
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Figure 6: log-log plots for |a2,4| versus a for various A in the slow quench regime. 
The fact that the plots tend to straight hues for positive values of log a means that 
02,4 scales as a power law for large a. Clockwise from the top left, the plots are for 
A = 7/3, 8/3, 11/3 and 10/3. 

9.3 Response for slow quenches 

In figure 5, we plotted 0(2A-4) as a function of r for the various values of A and a. The 
time is scaled by so that the different plots would equilibrate in approximately 
the same distance on the horizontal axis. As a grows large, the curves approach an 
inverted tanh graph, which is the expected adiabatic limit, i.e., 0(2A-4)(oo) 0(o)(t) — 
see eq. (8.2). Note that in this case there is no need to rescale the vertical axis since 
0(2A-4) is of the same order in all cases. 

As discussed in section 8, 02,4(00), which controls the entropy production, goes 
to zero in the adiabatic limit. Further, the analysis there showed that the leading 
contribution gave 02,4(00) oc 1/a for slow quenches — see eq. (8.10). This behaviour 
is revealed in our numerical results in figure 6. There log(— 02,4(00)) is shown as a 
function of logo; and we see that for large a, the results can be fit with a straight line 
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Figure 7: (Colour online) Plots of the deviation from the adiabatic response as a 
function of 0(o) for slow quenches with different speeds and A = 11/3 — see eq. (9.8) 
for the definition of 0(2A-4)- The curves correspond to a = 1 (grey), 2 (brown), 4 
(blue), 8 (purple), 16 (green), 32 (orange) and 64 (red). The dashed curve corresponds 

to &{O)0(O)- 



with a slope of approximately —1 in all the plots shown. Similar to the fast quench 
case, the straight lines are fit through the last three data-points in each plot.^ Further 
the intercepts of the straight lines in figure 6 should correspond to (minus the logarithm 
of) the coefficients given in eq. (8.11). We defer the detailed comparison of the results 
derived in section 8 and with the numerical simulations here until section 10. 

Again, although not shown, the same behaviour was also found for slow quenches 
in the case A = 4. This behaviour was also found to hold for A = 2 and 3 in [12]. 

For slow quenches, let us define the deviation of the response from the adiabatic 
limit (8.2) as 

4>{2A-4)ir) = a ( 0(2A-4)(1") - 0(2A-4) (oo) 0(0) (t) ) . (9.8) 

As discussed in section 8, this function should be approximately given by &(o)0(o); 
where 6(o) was the coefficient of the normalizable mode in the radial profile of the 1/a 
contribution. Figure 7 shows the deviation 0(2A-4) as a function of 0(o) for A = 11/3 
and different values of a. The dashed curve shows 6(o)0(o), where 6(o) was determined 
by the shooting method in section 8. As we can see in the figure, as a grows large, 

^Note that in the case A = 7/3 we fit the fine through three intermediate points, since our result 
for 02,4(00) contained a significant numerical error for the largest value of a shown. 
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Figure 8: (Colour online) Plot of 5 as a function of r for a = 1 and A = 8/3. The 
excitation time r^x and relaxation time r^q are shown as the first and final times, 
respectively, at which 5 crosses the threshold oi e = .05 (shown as the orange line). 

the deviation determined by our numerical simulations is converging on the expected 
curve. Those curves corresponding to larger a fit the dashed curve best. The curve for 
a = 64 lies practically on top of the limiting dashed curve. 

9.4 Excitation and relaxation times 

Next, we consider the excitation and relaxation times following the approach presented 
in [12]. First, we define 



i.e., the absolute value of the relative deviation of the response 0(2A-4) from the adia- 
batic limit 0(2A-4)(oo) 0(o). Then we define the excitation time t^x as the first time at 
which 5 reaches e = .05, where the latter was chosen as an arbitrary threshold. Simi- 
larly, the relaxation or equilibration time r^g is the latest time at which 6 drops below 
the e = .05 threshold. As an example, 6 is shown in figure 8 for a = 1 and A = 8/3. 
The vertical grid lines indicate the excitation time Tex and relaxation time Tgq. Note 
that our definitions of and Tgq is only expected to be meaningful for relatively small 
a. As a grows large, the response approaches the adiabatic profile (8.2) and so for 
sufficiently large a, S will never exceed the chosen threshold. 



S = 



0(2A-4) - 0(2A-4)(OO)0(O) 
0(2A-4)(OO) 



(9.9) 
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Figure 9: Plots of as a function of log a. The straight lines indicate the asymptotic 
behaviour for fast quenches. Clockwise from the top left, the plots are for A = 7/3, 
8/3, 11/3 and 10/3. 

Figure 9 shows the rescaled excitation time Ireil/o; as a function of log a for different 
values of A. For fast quenches {i.e., log a < 0), we see that \Tex\/oi approaches a 
straight line, indicating that the excitation time scales as /3a ct log a. Once again the 
straight lines are fit through the three points with the points corresponding to the 
fastest quenches, for each A. The constants can be determined as the slope of the 
fitted line in these plots. The plots also show that, as expected, the behaviour becomes 
irregular for logo; > 0, in particular for A = 10/3 and 11/3. 

Figure 10 shows the slopes of the straight-line fits in the previous plots as a function 
of A. This plot also includes data-points for A = 2 and 3 using the results in [12].^ 
The fact that = for A = 2 means that in this case t^^ has no log a-dependence 

^For these two points, we have used the expressions for \Tex\/(y. after the log (—log a) terms have 
been subtracted. 
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Figure 10: A plot of (minus) the slope of the fitted straight lines in figure 9. The 
datapoints lie approximately on the predicted line = A — 2 shown. Also 

shown are the data for A = 2 and 3 from [12]. 

for this type of quench. We find that the datapoints lie approximately on the line 

d\Tex\/0i 



9 log a 



= A-2, (9.10) 



which is numerically almost identical to the fitted line through all the data-points 
shown in the figure, namely 

5ke.|/« 1.003A-2.02. (9.11) 



9 log a 

Hence for fast quenches, i.e., a <^ 1, the excitation time scales as 

Te^~ (A-2)aloga. (9.12) 

Using the same approach, we also studied the relaxation time Tgg. In this case for 
all of the various A, we found 

Te, RiQ;°, (9.13) 

for fast quenches. That is, Tg^ is constant for small «, rather than scaling with a in 
some way, for all (fractional) A. This behaviour is consistent with the results for A = 2 
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and 3, found in [12]. In terms of the dimensionful time coordinate t, the relaxation 
time is t^g ~ l//i. In terms of the boundary theory then, the relaxation time is set by 
the thermal timescale for fast quenches. 

9.5 Behaviour of the energy and pressure 

Recall that our analysis was restricted to considering conformal dimensions in the range 
2 < A < 4. With this restriction, the change in the temperature AT and energy density 
AS will always be positive in our holographic quenches, as is evident from eqs. (6.25) 
and (6.26). On the other hand, eqs. (6.35) and (6.36) show that AT and AS may 
have either sign for the reverse quenches considered in section 6.3. In particular, in 
the adiabatic limit, 02,4(00) vanishes and so we have both AT < and AS < 0. Then 
as a becomes smaller, 02,4(00) grows and eventually AT and AS become positive. 
Specifically, eqs. (6.35) and (6.36) indicate: 

AT > for |a2,4(oo)| > ^ (A - 2f |0(2A-4)(-oo)| , (9.14) 

A£:>0 for |a2,4(oo)| > ^(A-2)|0(2A-4)(-oo)|. (9.15) 

Recall that 0(2A-4)(— c>o) = 0(2A-4)(oo) corresponds to the equilibrium response given 
in eq. (6.13). With our numerical simulations, we determined the value of a at which 
these thresholds are reached for various values of A and the results are shown in table 
2. We have also included the analogous results for A = 2 and 3 from [12]. Note the 
qualitative trend is that the threshold value of a grows monotonically as A increases. 
Note that eqs. (9.14) and (9.15) imply that the threshold for positive AS is greater 
than that for positive AT {i.e., larger a) for A > 2.5, while the thresholds are reversed 
for A < 2.5. Clearly, this behaviour is reflected in the results shown in table 2. 

Turning now to the change in the pressure AV, we have eqs. (6.27) and (6.37) 
for forward and reverse quenches, respectively. In this case, AV > in all forward 
quenches as long as A > 7/2 and in all reverse quenches for A < 7/2. Otherwise, 
the sign of AV will depend on the rate of the quench. Here we focus on the forward 
quenches with A < 7/2. In this case, AV < for slow quenches and as a decreases, 
the change in the pressure reverses its sign when 

AT' > for [02,4(00)1 > ^ (7 - 2A) (A - 2) [(/)(2a-4)(oo)[ . (9.16) 

The value of a at which these thresholds is reached for various values of A is shown 
in table 2. Again, we have also included an analogous result for A = 3 [12]. The 
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Table 2: Approximate upper bounds on a for which AT > and A£^ > for reverse 
quenches. The upper bounds on a for AT' > is for forward quenches when A < 3.5, 
with different values of A. For A = 2 the values of t^E and AT' are renormalization 
scheme dependent, see [12]. 



A 


AT 




AT' 


2 


0.58 






7/3 


0.68 


0.50 


0.23 


8/3 


0.77 


0.92 


0.68 


3 


0.86 


1.32 


1.32 


10/3 


1.00 


2.00 


5.5 


11/3 


1.41 


4.00 





qualitative trend is that the threshold value of a grows monotonically as A increases. 
Further, we may compare the above threshold to those in eqs. (9.14) and (9.15).^ In 
particular, the threshold for positive AP is greater than that for positive A£ for A > 3, 
while the thresholds are reversed for A < 3. We also find the threshold for positive AT' 
is less than that for positive AT when A < 2.75. Clearly, this behaviour is reflected in 
the results shown in table 2. 

10 Discussion 

In this paper, we continued the program initiated in [12] of studying quantum quenches 
in strongly coupled quantum field theories using holography. The process studied here 
was the quench of a four-dimensional conformal field theory made with a rapid tran- 
sition in the coupling of a relevant operator (9a from zero to some finite value A/. 
Ref. [12] had considered the special cases where the conformal dimension of the oper- 
ator was A = 2 and 3. Our holographic analysis allowed for operators with general 
conformal dimensions in the range 2 < A < 4. Through the gauge/gravity correspon- 
dence [5] , this quench was translated to a classical problem in five-dimensional Einstein 
gravity coupled to a negative cosmological constant and a (free) massive scalar field. 

''Note that here we are comparing a threshold for the forward quenches to thresholds in the reverse 
quenches. 
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In particular, the quench was implemented by introducing a time- dependent boundary 
condition on the scalar field in the asymptotically AdSs spacetime. Our discussion also 
only considered quenches of a thermal plasma with an initial temperature Tj in the 
boundary theory and was limited to a high temperature regime^ where A/ <C T^~^. 
Hence our calculations were perturbative in the ratio Xf/T^~^, which in the gravita- 
tional description meant that they were perturbative in the (dimensionless) amplitude 
of the bulk scalar. However, there was no restriction on the timescale At governing the 
transition rate of the coupling. In particular, with the transition profiles in eq. (7.1), 
our results were described in terms of the dimensionless parameter: a = TrT^At, i.e., 
the ratio of the transition timescale to the relaxation timescale of the thermal plasma. 
In our analysis, we paid special attention to the limit of adiabatic transitions with 
a — >■ oo and of very fast quenches with a ^ 0. A detailed discussion of the results was 
given in section 9. 

In all of these quenches, our computations implicitly gave the response in the one- 
point correlators of the relevant operator (Oa) and the stress tensor (Tij), as described 
by eqs. (5.4)-(5.6), to leading order in Xf/T^^^. In section 9, our discussion of the 
results was presented in terms of two gravitational parameters, (/)2A-4{t) and 02,4(00). 
With cq. (5.6), we see the first is directly related to (Oa) during the quench. We can 
rewrite this expression as 

using eq. (2.5) and various results from section 6. Further 02,4(00) controls the entropy 
production (6.24), as well as the changes in the temperature, energy and pressure as 
given in eqs. (6.25)-(6.27). One confirmation of our numerical simulations was that we 
found 02,4(00) < for all of our quenches with any values of a and A. With eq. (6.24), 
the latter ensures that the entropy production was always positive, in accord with the 
second law of thermodynamics. 

Another confirmation comes from our analysis of slow quenches in section 8. In this 
case with a ^ 1, it was shown that the linearized equation (3.5) for the bulk scalar 
can be solved using a power series in 1/a. While 02,4(00) vanishes for the leading 
adiabatic solution, we showed that a 1/a contribution appears at the next order in 
eq. (8.10). Using a shooting method, we explicitly solved for this contribution and the 
results were given in eq. (8.11). This same leading order contribution to 02,4(00) for 



(0.)-^(A-2)5^--^. (10.1) 



^In fact, this inequality is satisfied by the couphng and temperature throughout the quench. 
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slow quenches is addressed with the numerical results shown in figure 6. In these plots, 
the same approximate 1/a scaling was found with an asymptotic straight- line fit in 
logo; for a':^ 1, i.e., 

log 1 02,4(00)1 = — c — logo; . (10.2) 

The case A = 7/3 had the worst fit, with a slope differing from —1 by about 18%. 
We believe that the slow convergence in our numerical simulations for this case led to 
this relatively large error. The other three cases in table 3 had slopes that differed 
from the expected slope by 2% or less. As well as obtaining a fair match in the slope 
above, we can compare the intercept c coming from these numerical results with that 
calculated from the independently derived results in eq. (8.11). We see in table 3 that 
the intercepts derived from the two approaches agree very well. Recall that in figure 
7, we also showed that the full time-dependent profile of the numerical response (after 
subtracting the adiabatic profile) matched well with the form derived in section 8 for 
q; > 1. 

Table 3: Intercept in eq. (10.2) evaluated by two different methods: Cghoot is derived 
from the results in eq. (8.11), while Cnumer comes from fitting the data in figure 6. 



A 




^ II U IlK !' 


('.■^h nol ^'ri 'ti mer 
'■\^hool 


7/3 


3.93 


3.71 


5.6% 


8/3 


2.96 


2.98 


-0.68% 


10/3 


2.31 


2.26 


2.2% 


11/3 


2.22 


2.21 


0.45% 



In eq. (9.12), wc found an interesting scaling behaviour for the excitation time in 
fast quenches, namely r^x — (A — 2) a log a. On the one hand, this indicates that 
the excitation time is longer when the conformal dimension of the operator is larger 
but it also shows that Tex becomes shorter for faster quenches. In particular, Tex 
for a ^ 0, for which the quench profile (7.1) becomes a step-function at r = 0. In 
contrast, as shown in eq. (9.13), the relaxation time remains constant for a <^ 1. 
Hence independent of the precise values of A and a, the boundary system relaxes on 
the thermal timescale 1/T for fast quenches. 

Perhaps, the most interesting result coming from our analysis was the universal 
behaviour found in |02A-4| and 02,4(00) for a <^ 1. Of course, in terms of the boundary 
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theory, these results translate into universal behaviour in the response of (Oa) and in 
the thermodynamic quantities for fast quenches. First, the results in figure 1 indicate 
that the maximum value of |(62A-4| scales as a~^'^^~'^\ which with eq. (10.1), translates 
into a scaling for the expectation value of the quenched operator, i.e., 

max(e»A) cx (10.3) 

Beyond this scaling, figure 2 also indicates that the response and hence (Oa) approach 
a relatively simple universal form in the limit a — ?■ 0. These results are in agreement 
with those found previously in [12]. However, we might comment that our analysis here 
assumed that A was a fraction, whereas [12] studied the special cases A = 2 and 3. In 
both of these cases, the response also exhibited an additional contribution which scaled 
faster than shown in eq. (10.3) by an extra logarithmic factor. An extra logarithmic 
factor is also present for A = 4 [31]. 

The above scaling of |02A-4| also leads to the scaling of 02,4(00) oc a;~*^^^~^^ for fast 
quenches, as shown with the numerical data in figure 4. This contribution then domi- 
nates in eqs. (6.24)-(6.27) and so the changes of the various thermodynamic quantities 
induced by fast quenches exhibit the same scaling, e.g., 

for a 1. Again these results are in agreement with those found in [12]. However, the 
results there and in [31] indicate that eq. (10.4) is further enhanced by a logarithmic 
scaling for A = 2. While quenches by operators with conformal dimensions in the 
range 2 < A < 4 are covered by the analysis here and in [12], it would be interesting 
to understand if this universal behaviour extends to the allowed regime 1 < A < 2. 
Wc will consider this possibility in a later paper [31]. 

As noted in [12], given the scaling in eqs. (10.3) and (10.4), it appears that 'infinitely 
fast' quenches seem to be ill-dcfincd because physical quantities arc diverging as a — )> 0. 
Recall that in this limit, the quench profile (7.1) becomes a step-function at r = 0. 
Hence this issue is particularly notable since it is precisely such 'infinitely fast' quenches 
are studied in the seminal work on this topic [3]. However, we must contrast their 
description of a quench with the present approach. In [3], the system is evolved from 
t = —00 to 0~ to prepare the system in a far- from- equilibrium state of the 'quenched' 
Hamiltonian, i.e., , the ground state of the initial Hamiltonian. This state is then used 
as the initial condition at t = 0^ and the subsequent evolution of the system with the 
'quenched' Hamiltonian is studied. 
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Of course, since the present calculations are only perturbative in \f/T*~^, one can 
not take the singularities appearing in eqs. (10.3) and (10.4) for a — )■ too seriously. 
Hence it would be interesting to study the fast quenches by evolving the full nonlinear 
equations of the dual gravity theory. At present, our preliminary analysis suggests 
that in fact these singularities are physical [31]. In any event, the present holographic 
calculations illustrate that the gauge/gravity correspondence provides a versatile new 
framework for the study of quantum quenches. Undoubtedly, interesting new lessons 
will come from applying holography to study more general physical quantities and 
the behaviour of more complicated systems under a quench. This will help build our 
intuition for the behaviour of fast-changing quantum fields that occur when the external 
parameters are changed in laboratory experiments. 
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A Coefficients in the metric solution 

Here we list the expressions of the coefficients in the metric functions (3.11) and (3.12) 
in terms of the normalizable and non-normalizable modes of the scalar, 0(o) and 0(2A-4)- 
We only list the coefficients that are needed (and the ffist subleading coefficient) in 
calculating the boundary stress tensor and the expectation value of the operator (9a- 
Because 0,2,4 depends on it, we will give the expression for 02,5, even though it is 
subleading and has only vanishing contributions to physical quantities: 

02,5 = ^ (a(2A - 5)0(o)0(2A-2) - (2A - 3) (4 - A)0(o)</)(2A-2)) • (A.l) 
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The coefficients of tlie terms witfi negative powers of 2 A in Ap are given by 



(4 - A)(/)2 

^'^ 3(7 -2A) ' ^ ^ 

(2A2 - 6A + 15)0fo) + (2A2 - ISA + 24)) 0(0)0(0) 
"^'^ - 12(A-3)(9-2A) ' ^^-^^ 

(4 - A) (-3(A - 2)(9 - 2A)0(o)0(o) + (2A2 - 15A + 36))</.(o) 0(0)) 
"^''^ " 36(5- A)(A-3)(9-2A) 

a2,8 = -(l2(5-A)(4-A)2(A-3)2(2A'-llA + 30)0^0) 
-3(7 - 2A)2(5 - A)(2A3 - IQA^ + 56A - 54)0jo) 
-4(A - 3) (7 - 2A)(2A^ - 14A + 21) (2A^ - 17A + 39)0(o)0[o) 
+(4 - A) (A - 3) (9 - 2A)(2A - 7)(2A2 - 17A + 51)0(o)0|o5) 
/ (288(5 - A) (4 - A) (A - Sfi^A'' - 36A + 77)) . (A.6) 

Of tfie coefficients to tfie terms witli positive of powers 2 A in Ap, the coefficient 

^^'^ - 6(2A - 1) ^^-^^ 

will be subleading. 

The coefficients in Ep corresponding to integer powers of p are given by 

^ _A(4-A)0(o)0(,^_4) 

2,5 36 ' ^ ' 

S2,6 = (A(5- A)0(o)0(2A-4) + (4-A)(A+l)0(o)0(2A-4)) • (A.9) 
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The coefficients of tlie terms witfi negative powers of 2A in Ep are given by 



^^'^ ^ -12(7-2A)' ^^-^^^ 
^ _(5-A)0(o)0(o) 
^'^ 6(9 -2A) ' ^ ^ 

^ 2(5 - A)^(A - 3)0^,) + (6 - A)(4 - A)(7 - 2A)0(o)0(o) 
"^''^ 24(5- A)(A-3)(9-2A) ' ^ ^ ^ 

3(6-A)(5-A)(7-2A)0(o)0(o) 

72(5- A)(A-3)(11 -2A) 
, (7-A)(4-A)(9-2A)0(o)0gj _ 



72(5- A)(A-3)(11-2A) 
(72,9 = -(l2(8 -A)(4-A)2(A- 3)202^)- 3(7 -2A)2(6-A)202^) 

-8(7 - A) (5 - A) (A - 3) (9 - 2A)0(o)0g; 

-(8 - A)(A - 3)(11 - 2A)(9 - 2A)0(o)0g;) 

/(^576(6 - A)(A- 3)2(11 -2A)). (A.14) 
Of tlie coefficients to tlie terms witfi positive of powers 2 A in Ep, tfie coefficient 

^2 5 = ^^1^ (A.15) 
12(2A - 1) ^ ^ 

is subleading. 

B Coefficients in the Fefferman- Graham coordinates 

B.l The time and radial coordinates 

Tfie expansion of tfie EF time and radial coordinates in terms of the FG time and radial 
coordinates was given in eqs. (4.5) and (4.6). The expressions of the coefficients can be 
written in terms of the coefficients 0(o), 0(2A-4) and 02,4. First, the terms containing 
Vn and pn are given by the series 

„ 3rV ll^V 23riV^ ^r. ■,^ 

> = -r h . . . , (B.l) 



n 



128 1024 



48 



Next, the terms of order £^ are given by 

1 5 5 1 

r ^ "^'^ + 10 



n=5 

= ^^^^"^"2,4 + Tf^rV^ (/xa2,5 - 5*02,4) + ... , (B.4) 

n=5 



where the terms of order only make subleading contributions in the calculated quan- 
tities. The coefficients of the terms with factors of r"^^ in r//x are given by 

(7 - 2A)fj,^^-^a2,^ 
~ 4(4- A)(9-2A) ' 

^8-2A ^4^2 _ + 55) dta2A + 4(4 - A) ^^2,5) 



i^2 = A* 

1/3 = -II 



8(5- A)(ll- A)(9-2A) 
_2A (4A2 - 32A + 61) d^a2,4 
8(5-A)(ll-2A)(9-2A) 
-2 (4A2 - 34A + 71) ^^2,5 + 2(9 - 2A) ^^^2,6 
8(5-A)(ll-2A)(9-2A) ' 
8_2A (4A2 - 34A + 67) d^a2,4 - 6(2 A^ - 18 A + 39)iidfa2,5 



48(6- A)(5-A)(11-2A) 



_ 8_2A 6 (4A2 - 38A + 89) pL''dta2fi - 24(5 - A) ^=^0^2,7 . 

^ 48(6- A)(5-A)(11-2A) ' ^ ' 

where a2,n are the coefficients defined in eq, (3.11) and given explicitly in eqs. (A.2)- 

(A.5) — implicitly, functions of the FG time t here. Similarly, the coefficients in p with 
factors of r"'^^ are given by 

- 2(9 -2A) ' ^^-^^ 
(5|«2,4 - 2iidta2,5 + 2ij?a2fi) 

^ ' ^ ^ ^ 

12(11 - 2A) ' ^^-^^ 



e 



4 



^9-2A ((2A2 - 13A + 30) ii^a2,4 - (4 - A)a| 0^2,4) 
32(6- A) (4 -A) 
^4^dfa2,5 - I2fi^d^a2,(i + 24//3a,a,^, - 24//^a2,8) 
96(6 - A) 



(B.IO) 



49 



The coefficients of tlie terms witli factors of r"^^ in v are given by 

- 4A(2A + 1) ' ^^-^^^ 

- 8A(A + 1)(2A + 1) ' ^^-^^^ 

2^ (4A2 - 3) 9^/32,4 

^3 — // 



8(4A3 + 12A2 + llA + 3) 
,^ 2 (4A^ + 2A - 1) - 2(2 A + 1)^^/^2,6 

^ 8(4A3 + 12A2 + llA + 3) ' ^ ^ ^ 

_ 2A (4A^ + 2A - 5)91^2,4 - 6(2A^ + 2A - l)fxdfP2,5 
^ 48(A + l)(A + 2)(2A + 3) 

(4A^ + 6A + 1) /.^a,^2,6 - 4( A + 1)^3^2,7 
^ 12(A + l)(A + 2)(2A + 3) ' ^ ' 

wliere ^2,n are the coefficients defined in eq. (3.11) — implicitly, functions of the FG 
time t here. Similarly, the coefficients in p with factors of r^^ are given by 

Ci = ^4/^; (B.15) 

^ 2(2ATT^ ' ^^-^^^ 

^ 8(ATI^ ' ^^-^^^ 

^ 12(2A + 3) ' ^^-^^^ 

_ ,^^, 3 (2A^ - 3A + 10) /x^/32,4 - Aat'/g2,4 
^ 96A(A + 2) 

2A+1 A {ipdf(32,, - 12/^^9,^^2,6 + 24/x=^a,^2,7 - 24/X^^2,8) . 
96A(A + 2) ■ ^^-^^^ 

B.2 The metric 

The expansion of the metric in the FG coordinates was given in eq. (2.8). The nonzero 
components for our purposes are the parts Gqq, which correspond to the EF metric 
function —A, and the diagonal components of G^, which correspond to E^. The order 
£^ terms in the metric are given by 

glS + r'giS - -l + ^ + o(r«), (B.20) 

a^+r'glP = l + 4^ + o(r«), (B.21) 
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where in the second hne, we are indicating the three individual diagonal components, 
i.e., there is no implicit sum over i. Next we list the terms of order 

3 

C(4)oo = -^/^%,4, (B.22) 

C(4)ii = (^a2,4 + ^A(4- A)0(o)0(2A-4)^ • (B.23) 

The coefficients of the terms with factors of r"^^ in Gqo are given by 

_ (7-2A)/-^V, . 

«(o)oo 2(4 - A) ' ^ 

J o,.8-2a (A - 3)9^0:2,4 + (4 - A)//Q;2,5 . , . 

0(1)00 = -2/i — —r , (B.25j 



9 - 2A 

•^(2)00 = 1^ 



s_,^ (A-3)(7-2A)a,^«2,4 
4(5- A)(4- A) 
2a2(7 - 2A)fidta2,5 - 2(9 - 2A)fi''a2,6 



+^ 4(5^A^ ' ^^-^^^ 



C^(3)00 = -/^^ 



(A - 3) (7 - 2A)gfQ2,4 + 3(4 - A) (7 - 2A)/^a|a2,5 
3(4A2 - 40A + 99) 



«-2A (4 - A)a,a2,6 - 2(9 - 2A)/x^ ((5 - A)/.a2,7) . ,^ 
+^ 4A2-40A + 99 ' ^^-^^^ 



t(4)00 — —A* 



-2A 



(2A3 - 25A2 + 107A - 162) fjL^a2,4 



8(6- A) (5 -A) 
.,^ 4(4 - A)(7 - 2A)^dfa2,, - 12(9 - 2A)(4 - A)fi''d^a2,e 

48(6- A)(5- A) 

8-2A 



24(9 - 2A)n^dta2j - 24(11 - 2A)/x^a2,8 



48(6- A) 



.2a(A~3)(7-2A)9,V4 



I..8-2A V" "7V' -"7'^f^:^,4 m OS'l 

+^ 48(6-A)(5-A) ' ^^-^^^ 
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where a2^n are the coefficients defined in eq. (3.11) — imphcitly, functions of the FG 
time t here. Similarly, the coefficients in Ga with factors of r~^^ are given by 



d{o)n - -/^ 2(4 - A) ' (l5-^yj 

, 8-2A ^t«2,4 - 2(9 - 2A)dta2,5 - fia2,5 + 2(9 - 2A)/x(72,6 on^ 

"(l)n = Q _ 2A ' (D.dUj 

> „8-2A^I«2,4-4(5- A)92(T2,5 

"(2)ii = -/^ 



4(5 - A) 

A)/xa 

2(5- A) 



8-2aM«2,5-4(5- A)/i5ta2,6-/i^a2,6 + 4(5-A)/xV2,7 



^^(3)ii = _ (^'«2,4 - 2(11 - 2A)a3a2,5 - Si,d^a2,5 + 6(11 - 2A)i,d^a2,e 

+9*02,6 - 2(11 - 2A)dta2,7 - O/x^ (^^2,7 + 2(11 - 2A)/i(T2,8) ) , (B.32) 

where (T2,n are the coefficients defined in (3.12) — implicitly, functions of the FG time 
t here. 

The coefficients of the terms with factors of r^^ in Gqo are given by 

(2A-l)/.^^/32,4 . 
e(o)oo = ^ ; (B.33) 

_ 2a (A- W2,4-A/.^2,5 .^,4. 

6(1)00 - 2/i 2A+ 1 ' ^ 

where ^2,n is defined in eq. (3.11). Similarly, in Gu the coefficients of the r^^ are given 

by 

- ..2a /^2,4 - 4A6'2,5 . orN 

e(o)n - -/^ , {B.Sb) 

= ^ 2A+ 1 ~' ^ ^ ^ 

where T2,n was defined in eq. (3.12). 

B.3 The scalar field 

The coefficients of the scalar field /(„) and in the FG coordinates can be written 
in terms of the coefficients 0(o) and 0(2A-4) in eq. (3.9). The coefficients of r^~^ in $p 
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are given by 

r2(7-2A)9f0(o) r=^(9 - 2A)af 0(o) 



4(A-3) 12(A-3) 
,fl^^ . (4A2-40A + 99)9,V(o)\ , 

(B.37) 

while the only coefficient of that may play a role is 

5^(0) = //'^0(2A-4). (B.38) 

C Boundary stress-energy tensor and {Oa) 

In section 5, we constructed the holographic action Sreg — Sf,uik + Sgh by + S count which 
is now finite, i.e., no divergences appear in the limit e — > 0. Hence, so are all quantities 
that can be calculated from this action. Of course, the latter includes the one-point 
functions of the stress tensor and operator Oa- In order to calculate these expectation 
values, we need to vary Sreg with respect to the boundary metric and the coupling to 
boundary operator, respectively. Recall that the boundary metric g^^^ appears as the 
leading coefficient in the expansion (4.7) of the bulk metric. Similarly, the coupling A is 
proportional to the leading coefficient /(o) in the expansion (4.8) of the bulk scalar. We 
will establish our conventions for the precise normalization of the coupling in section 
6.1 and so at this point, we simply introduce a dimensionless proportionality constant^ 
with £/(o) = ax A. Then the desired one-point functions are given by [24] 



lim 

e->0 














lim 




1 ^Sreg 


Kb 



5Sreg 5^ 



(C.l) 



and 



lim 








5X 


lim 




"A f ' 







5^ab 5f{0) 5^ Sf^o) 



+ ^--^ . (C.2) 



^The constant a\ is fixed in eq. (6.16). 
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We proceed by first evaluating the variations of the action with respect to 7 and 
then the variations of 7 and $ with respect to g^j^^ and /(o). The variation of the 
action on the cut-off surface r = e is: 

V-7 ^^nh 2 



'-7 



^ ('S'bulk + 'S'ghBy) 



'-7 



count __^a'> 

(57a6 2 



6 + 



4- A 



1 



IGTrCjy'' (55'count 



-7 



5$ 



24 (A - 3) 
= - (4 - A) $ - 



4 (A - 3) 
(V" V'' - 7"^ □) - 67"^7'"^9c$ 

1 



2 (A - 3) 



(C.3) 



In the above we only showed the terms that would contribute at the orders of £ we are 
considering. Using the results of appendix B, we find the variation of the bulk fields 
with respect to their boundary values (at leading order) are 



X (0) 
hah 



5^ 

Kb 



-g^a6^V(0)r' 



6-2A 



Ir 



4-A 



(C.4) 



The second equation in eq. (C.4) leads to a vanishing contribution to the stress tensor 
in eq. (C.l). The third equation above only contributes to eq. (C.2)when A = 4 
and hence can be ignored for our purposes. Inserting the variations from eqs. (C.3) 
and (C.4) into eqs. (C.l) and (C.2), we get the results presented in the main text in 
eqs. (5.4)-(5.6). Note that for these final results, we have replaced /(o) and 5f(o) by their 
expressions in terms of 0(o) and 0(2A-4) given in appendix B.3. 
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